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The SUSY signals in the dominant stau-neutralino coannihilation region at a 500(800) GeV Hnear 
collider are investigated. The region is consistent with the WMAP measurement of the cold dark 
matter relic density as well as all other current experimental bounds within the mSUGRA framework. 
The signals are characterized by an existence of very low-energy tau leptons in the final state due to 
small mass difference between fi and x\ (5-15 GeV). We study the accuracy of the mass difference 
measurement with a 1° active mask to reduce a huge SM two-photon background. 



1 Introduction 

The recent measurement of cold dark mat- 
ter (CDM) relic density from WMAP^ along 
with the Higgs mass bound and the b — > 
S7 constraint have restricted the parameter 
space significantly^ within the framework of 
minimal supergravity (mSUGRA) model. '^'^ 
One prominent parameter space is the region 
where the mass difference (AAf) between the 
lighter stau (fi) and the lightest neutralino 
is about 5-15 GeV. This small mass dif- 
ference allowed the fi to coannihilate in the 
early universe along with the x\ in order 
to produce the current amount of the CDM 
(Xi)- The coannihilation region has a large 
extension for mi/2 up to 1-1.5 TeV, and can 
be explored at the LHC. The main difficulty, 
however, in probing this region is to detect 
very low-energy taus in the final state of the 
SUSY events due to the small AAf value. 

In this paper, we report a feasibility 
study of measuring the small mass difference 
in this Ti-x\ coannihilation region at a 500 
GeV linear collider (LC). 

2 mSUGRA Parameter Space 

The mSUGRA model depends on only four 
parameters and one sign. These are mo (the 
universal soft breaking mass at the GUT scale 
Mq); nil 12 (tbc universal gaugino soft break- 
ing mass at Mq)\ Aq (the universal cubic soft 
breaking mass at Mq); tan/3 = {H2) / {Hi) at 
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Figure 1. Allowed region in the mo-mi/2 plane from 
the relic density constraint for tan /3 = 40, Aq = 
and n> Q. The details are provided in text. 



the electroweak scale; and the sign of /i, the 
Higgs mixing parameter in the superpotential 
{W^=tiHiH2). 

Figure 1 is an example of the allowed re- 
gion in the mo-mi/2 plane for tan/? = 40 
with j4o = and > 0. The most impor- 
tant experimental results for limiting the pa- 
rameter space are: (1) The light Higgs mass 
bound (pink lines in the figure) from LEP:^ 
Mh > 114 GeV; (2) The b ^ s -I- 7 branching 
ratio (brick red region):® 2 x 10^** < B{B — > 
A^7) < 4.5 X 10--*; (3) Previous CDM (x?) 



< 0.21 (yel- 
low band) from balloon flights (Boomerang, 



density bounds of 0.07 < fl~oh'^ 

Xi 



Maxima, Dasi, etc.) and the 2(t bound of 
0.095 < n~oh^ < 0.129 (blue band) from 

WMApi; (4) The bound on the hghtest 
chargino mass:^ xf > 104 GeV; (5) Possible 

muon magnetic moment anomaly (light blue 
region to be excluded if Sa^ > 11 x 10~^°).® 
It is striking to learn that only two SUSY 
production processes can be studied at a 500 
GeV LC: f^f^ and X2Xi- The kinematical 
reaches via the fj^fj" and X2X1 production 
are also shown in Fig. 1. The maximum reach 
in mi/2 along the coannihilation band can 
be expected via e+e" t^t^ — > (r+Xi) + 

We use the hadronic final state of tau 
(rh ) since it has larger branching ratios. Due 
to the small AM value, the taus in the final 
states are low energy and hence harder to de- 
tect. 

3 SUSY Signals at 500 GeV LC 

In order to optimize the event selection cuts, 
we choose three points of mo = 205, 210 and 
220 GeV for mi/2 = 360 GeV, tan/3 = 40, 
> 0, and Aq = 0. The SUSY masses 
given by ISAJET^ are summarized in Table 1. 
There are two major Standard Model (SM) 
background processes: (i) four-fermion final 
state i>iyT^T~ arising from processes such as 
diboson (M^M^, ZZ) production, and (ii) two- 
photon processes e+e~ 7*7* + e+e~ 
T^T~ (or qq) + e+e~ where the final state 
e+e~ pair are at a small angle to the beam 
pipe and the jets fake a r+r~ pair. 

The production cross-sections for SUSY 
(ISAJET) and SM four-fermion processes 
(WPHACTiO) are listed in Table 2 for a 500 GeV 
LC. We choose with right handed (RH) po- 
larized electron beams to enhance the fj^ff 
events over the X1X2 four-fermion 
events. 

In Table 3, we summarize the event se- 
lection criteria for the RH case. The; Monte 
Carlo (MC) events are generated, simulated 
and analyzed using the following programs: 



Table 1. Masses (in GeV) of SUSY particles in three 
representative scenarios of AM for mi/2 = 360 GeV, 
tan/3 = 40, /i > 0, and Aq = 0. 



MC mo 
Pt. 


M~o M~o AM 
X2 ^1 Xi 


1 205 

2 210 

3 220 


274.2 147.2 142.5 4.7 

274.2 152.0 142.5 9.5 

274.3 161.6 142.6 19.0 



Table 2. SUSY and SM production cross sections {cr ■ 
B(t — > Tji in fb) for polarization for electron beams 
of P(e-) = -0.9(RH). 



SUSY Pt. 1. 


X?X^ 


0.43 




^+^- 
'1 '1 


28.2.') 


SUSY Pt. 2. 


X?X^ 


0.39 




f+ff 


25.85 


SUSY Pt. 3. 


x?x^ 


0.38 




f+ff 


22.95 


SM (four fermion process) 


7.84 



ISAJET'^ to generate SUSY events; WPHACT1° 
for SM backgrounds; TAUDLA^-^ for tau de- 
cay; a LC detector simulation^ to reconstruct 
jets with JADE algorithm. In our calcula- 
tion, beamstrahlung and bremsstrahlung are 
included in both ISAJET and WPHACT. 

The accepted number of signal and back- 
ground events are summarized in Tables 4 
and 5. It should be noted that the number 
of SM 77 events with the forward electrons 
just below 3° are 11400. The acceptances 
for f+ff events are 11.2%, 5.9%, and 0.86% 
for AM= 19, 9.5, and 4.7 GeV, respectively 
with l°mask. The acceptance drops fast as 
AM goes below 5 GeV. For example, 0.23% 
for AM= 3.8 GeV (m,o = 204 GeV). We see 
the robust discovery significance for the sig- 
nal events for AM> 5 GeV with 1° mask in 
Table 6. We conchide that the mask is essen- 
tial to detect SUSY events in this region of 
parameter space. 
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Table 3. Event selection criteria for the RH {V = —0.9) case. 



Variable(s) 


Cuts 


Nj,t{Ejet > 3 GeV) 
Th ID 


2 

1, 3 tracks 


.Jet acceptance 


|cos((9,e/)| < 0.65 

-(.).() < COsO{j2.Pr,»)] < 0.(j 


Missing pt{^t) 


> 5 GeV 


Acoplanarity 


> 40° 


Veto on EM clusters 
or electrons 


No EM cluster in 5.8° < ^ < 28° with E>2 GeV 
No electrons within 6 > 28° with pT > 1.5 GeV 


Beam mask (l°(or 2°) - 5.8°) 


No EM cluster with E > 100 GeV 



Table 4. Number of SUSY events expected with 500 

fb-i 



Table 6. Significance {Ns/\/Nb) with 500 fb'^ for 



for the RH 


case. 








SUSY discovery using 1° mask. 






Process 


AM 


= 4.7 


9.5 


19 




Process (RH) 


AM^ 4.7 


9.5 


19 


A2A1 




15 


26 


29 




f+f- 


10 


63 


101 


f+f- 




122 


786 


1283 









Table 5. Number of SM events expected with 500 
fb-i. 



SM foTir-fermion 


129 


SM 77 2-5.8° Mask 
1-5.8° Mask 


248 
2 



4 Mecisurement of Stau Neutralino 
Mass Difference 

Since AM is small, it needs to be measured 
with a very good accuracy. We choose the 

invariant mass Mcff = M{ji, j2, ^ of two 
r-jets and missing energy as a key discrimi- 
nator. We generate high statistics MC sam- 
ples for the SM and various SUSY events (by 
changing the mo value) and prepare the tem- 
plates of the Meff distributions for the SM, 
XiX2> and f^f^ events. 

We then generate the MC samples equiv- 
alent to 500 fb"^ of luminosity for particular 
AM values and fit them with the template 
functions. For example, in Fig. 2 we show the 
fitting of the 500 fli^^ MC samples for Point 
2 with the templates for toq — 210 GeV and 



calculate the of the fits. Here the value 

is calculated as x = ^ ~ 1 

where Ni is the number of events in z-th MeflF 
bin of the 500 fb^^ sample, CjF/ is the cor- 
responding value for the template "j" where 
j is for SM, f+f]- or X1X2 processes. Cj is a 
normalization parameter and a free variable 
except for the SM process. This is because 
we should be able to measure the SM events 
very well before we discover SUSY events. 

We scan the range of mo = 203-220 GeV 
and plot the Ax^ = — Xmin Fig- 3- The 
Ax^ value is minimum for the template for 
mo = 210 GeV. We find that la in the Ax^ 
corresponds to 9.5 ± 1 GeV, where the true 
value of AM for the Point 2 is 9.53 GeV. 

We repeat the same study for different 
stau masses i.e. for different AM values and 
two different beam mask designs (1° and 2°). 
For AM 5 GeV, a beam mask of 1° is 
crucial. The accuracy of mass determination 
for is summarized in Table 7, showing the 
uncertainties are at a level of 10%. 
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fltting mQ=210 'data' with 
mn=210 templates 



X1X2 

Vi 

SM 

combined 




X^/n.d.f. 

N(Xi5C2) 

N(SM) 



70.18 / 98 
16.4 ± 16.6 
764.8 ± 30.6 
377.6 (fixed) 
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Figure 2. M^g (= M{ji, j2, f!)) distributions for a 
500 fb-i MC samples for SUSY (mo = 210 GeV) 
and SM events, being fitted to the templates for mo 
= 210 GeV. 
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Figure 3. (defined in the text) for Point 2 as a 
function of AM. 
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Table 7. Accuracy of the AAI determination for dif- 
ferent beam mask designs. "-" means we cannot de- 
termine with 500 fb~^. 



AM 


(500 fb"^) 


AAf("500 fb~^" expt.) 


2° mask 1° mask 


4.76 
9.53 
12.4 
14.3 


122 
787 
1027 
1138 


4.74t-^ 
9.5tl:l 9.5tl:0 



5 Conclusion 

At 500 GeV LC, it is crucial to instrument an 
active mask to detect very forward electrons 
down to 1° for measurement of the small 
AM. The expected accuracy is 10% (20% 
for AM - 5 GeV) with 500 fb"\ 
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